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Confinement  Properties  of  Non-neutral  and 
Neutral  Plasmas  in  an  Axially  Symmetric  System 


Teruo  Tamano  and  Isao  Katanuma 

( 'nNcrsity  ofTsukuba,  Tsttkuba,  Ibaraki  305-0045,  Japan 


Abstract.  A general  approach  is  established  to  examine  confinement  properties  of  both  non-nentral  and 
neutral  axially  symmetric  plasmas.  In  this  approach,  problems  are  reduced  to  finding  total  minimum  energy 
states  under  tire  conservation  of  total  angular  momentum  and  total  number  of  particles.  First,  the  well- 
known  result  of  a non-neutral  plasma  in  a cylindrical  system  is  reproduced  for  checking  the  approach. 
Second,  the  approach  is  applied  to  analysis  of  confinement  properties  of  tokamak  plasmas  with  a transport 
barrier  and  it  is  shown  that  the  tokamak  with  an  internal  transport  barrier  is  in  the  minimum  energy  state. 


INTRODUCTION 


It  is  well  known  that  a single  species  plasma  in  a Penning  trap  has  a very  long  lifetime. 
This  long  confinement  can  be  understood  by  the  constancy  of  the  total  canonical  angular 
momentum  in  a cylindrical  system  [1,2]  The  total  canonical  angular  momentum  is  given 
by  Po  = 2j(mju0fj  +e  jA)j’'j),  where  the  first  term  is  the  kinetic  term  of  the  angular 
j 

momentum  and  the  second  term  is  the  vector  potential  term  of  the  angular  momentum.  In 
the  case  of  a single  species  plasma,  the  vector  potential  term  is  much  larger  than  the  kinetic 
term.  Then,  the  total  angular  momentum  in  a uniform  magnetic  field  can  be  written  as 


J 


Accordingly,  the  location  of  the  single  species  plasma  particles  is 


constrained,  and  this  is  considered  as  the  essence  of  the  long  confinement 


On  the  other  hand,  in  a case  of  a neutral  plasma,  pairs  of  ions  and  electrons  can  escape 
together  from  the  system  without  changing  the  vector  potential  term  of  the  total  angular 
momentum.  This  is  considered  to  be  the  reason  why  the  confinement  of  neutral  plasmas  is 
poor  However,  the  situation  changes  when  the  mechanical  part  of  the  total  angular 
momentum  dominates  over  the  vector  potential  term.  Then,  the  constancy  of  the  total 
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canonical  angular  momentum  places  a constraint  similar  to  that  for  the  single  species 
plasma. 

In  this  paper,  we  establish  an  approach  to  examine  confinement  properties  of  both  non- 
neutral and  neutral  plasmas  under  the  constraint  of  conservation  of  the  total  angular 
momentum.  First,  the  approach  is  applied  to  a cylindrical  single  species  plasma  in  order  to 
check  its  validity.  Then,  the  approach  is  applied  to  a tokamak  plasma  with  a transport 
barrier. 


CYLINDRICAL  SINGLE  SPECIES  PLASMA 

In  the  case  of  a cylindrical  single  species  plasma,  the  total  angular  momentum  Pq  and 
the  total  number  of  particles  N per  unit  length  may  be  given  by 


Po 


2jt 


:Jo  n(mtl°  + eA° )? ‘ 


dr. 


(1) 


The  total  energy  U per  unit  length  can  also  be  conserved,  where  U is  the  sum  of  the  total 
kinetic  energy  K and  the  total  potential  energy  W,  i.e.,  U =K  + W The  total  kinetic 
energy  and  the  total  potential  energy  can  be  written  as 


K 


= nmuo  rdr , W = 


2 nf  (— nkT  + — e0E2  + — — )rdr . 

Jov2  2 ° 2/«() 


(2) 


The  plasma  equilibrium  under  the  conservation  of  the  total  angular  momentum  and  the 
total  number  of  particles  can  be  obtained  by  finding  the  state  for  the  minimum  total 
potential  energy  W . This  problem  is  equivalent  to  finding  the  state  for  the  maximum 
kinetic  energy  K when  the  total  energy  U is  conserved 
The  force  balance  equation  is  given  by 


E + ugB  =- 


1 d{nkT) 
en  dr 


(3) 


where  E is  the  radial  component  of  the  electric  field  Eliminating  ug  from  Eq.  (2)  and 
Eq  (3),  one  rewrites  the  total  kinetic  energy  as 
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ra  ,r  1 d(nkT)'  2 , 

K = ji  I nm(F.  + ~)  rdr . 

Jo  cn  dr 


(4) 


Since  the  vector  potential  term  of  the  total  angular  momentum  for  a single  species 
plasma  is  much  larger  than  the  kinetic  term,  the  total  angular  momentum  in  a uniform 
magnetic  field  is  given  by 


Pa  = neBr^dr. 


(5) 


By  use  of  the  Poison  equation 


n g»  1 d(rE) 
e r dr 


(6) 


Equation  (5)  becomes 


dr 


(7) 


Performing  a partial  integral  of  Eq.  (7),  one  obtains 


P0  = JTt[)It{E(a)a2  -2 f{)Er2dr). 


(8) 


The  total  kinetic  energy  Eq.  (4)  can  also  be  rewritten  by  use  of  Eq.  (6)  as 
f d(rE) ( r 1 d(nk'f)^ 


Tce()m  jid(rk) 


v_+ — -)  dr 

dr  en  dr 


(9) 


Performing  a partial  integral,  one  gets 


K .!x.\\fme  t jLf «2*n>f  _ fo-Et-oi 

eBl  |^[  en  dr  J()  *Ki  dr  cn  dr  j 


(10) 


Accordingly,  the  problem  is  reduced  to  finding  the  maximum  value  of  Eq  (10)  under  the 
constraint  of  conservation  of  the  total  angular  momentum  given  by  Eq.  (8),  This  variation 
problem  is  solved  by  using  the  Lagrangian  undetermined  multiplier  method: 
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5(K  -APe)-  0, 


(11) 


where  A is  an  adjustable  constant.  Equation  (11)  leads  to 


d_ 

dr 


{^eAe 

dr 


en  dr 


(12) 


Integrating  Eq.  (12)  and  using  the  condition  uq(G)  = 0,  one  obtains  the  following  relation 


„ 1 d(nkT) 

E+ — ^ 

en  dr 


Substituting  Eq.  (13)  into  Eq.  (6),  one  gets 


d r d(nkT) 
dr  en  dr 


(13) 


(14) 


A good  approximate  solution  for  the  plasma  density  is  given  by  n = n{  0)(  1 - e(r  n)j  '>'lJ  ) , 

where  the  Debye  length  ^■^eQkT(0)/e2n(0).  This  represents  the  well-known 
equilibrium  with  a density  profile  nearly  constant  up  to  the  plasma  radius  a and  sharply 
falling  down  to  null  within  a few  Debye  length  [1,2]. 


TOROIDALLY  SYMMETRIC  TOKAMAK  PLASMA 

In  the  case  of  a neutral  plasma  in  a toroidally  symmetric  system  such  as  a tokamak,  the 
total  toroidal  angular  momentum  Pq  is  conserved.  For  flux  surfaces  having  a concentric 
circular  cross-section,  one  may  write 

Pfj,  = 4jt2 Rq2 nmn^rdr , N = 4jz2  Ryj^imlr . (15) 

where  Aq  is  the  major  radius  of  the  axis  and  N is  the  total  number  of  particles.  The  total 
energy  U , the  sum  of  the  total  kinetic  energy  K and  the  total  potential  energy  W,  may 


559 


also  be  conserved:  U = K + W . The  total  kinetic  energy  averaged  over  flux  surfaces  may 
2 2 2 2 2 

be  written  as  K =4ji  J^f  nnni  n/r,  where  u - u#  +Uq  . The  total  potential  energy 
averaged  over  flux  surfaces  may  be  written  as 

2 

W = 4^2/^{2wA(7;  + 7;,)  +y0E2  + ~}rdr.  (16) 

Thus,  the  maximum  total  kinetic  energy  corresponds  to  the  minimum  total  potential  energy. 
Therefore,  the  problem  of  finding  the  minimum  potential  energy  state  is  reduced  to  find 
the  maximum  kinetic  energy  state  under  the  constraint  of  the  conserved  total  density  and 
total  angular  momentum  This  solution  is  again  obtained  by  the  Lagrangian  undetermined 
multiplier  method: 


8{K-X{N  -A27^)  = 0, 


(17) 


where  Aj  and  A2  are  constant  Equation  ( 17)  can  be  written  in  the  form 

(u2  - A!  -A2Mf)-J-+(2K*-A2)«-0.  (18) 

dtl/p 

The  solution  for  Eq.  ( 1 8)  is  given  by 


u - Aj  - A2  iifp 


(19) 


COMPARISON  OF  TEORETICAL  RESULTS  WITH  DIII-D  DATA 

The  theoretical  result  is  tested  against  actual  tokamak  experimental  data.  DIII-D 
plasmas  with  an  interna!  transport  barrier  (ITB)  are  chosen  as  a test  base.  Density  profiles 
of  the  DIII-D  plasmas  are  computed  from  the  measured  toroidal  angular  (rotation) 
velocities  based  on  Eq.  (19).  Figure  1 shows  the  measured  density  profile  of  QDB 
(Quiescent  Double  Barrier)  mode  DIII-D  plasma  as  the  solid  curve  [3]  together  with  two 
examples  of  calculated  theoretical  density  profiles  (dashed  and  dotted  curves).  Either  of 
them  are  a very  good  fit  except  for  some  mismatches  near  the  edge. 
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DIII-D  QDB  Density  Profile  Fit 


FIGURE  1.  The  normalized  density  profile  of  DIII-D  QDB  (shot  number  106919  at  3 s)  shown  by  the  solid 
curve  and  two  theoretical  fits  shown  by  dashed  and  dotted  curves. 


This  indicates  that  the  DIII-D  plasma  with  an  internal  transport  barrier  is  indeed  a 
minimum  (potential)  energy  state  under  the  constraint  of  the  conservation  of  total  angular 
momentum. 


CONCLUSION 

An  approach  is  established  to  examine  equilibrium  profiles  for  both  non-neutral  and 
neutral  plasmas  under  the  constraint  of  conservation  of  the  total  angular  momentum.  One 
important  result  obtained  from  this  examination  is  that  a tokamak  plasma  with  an  internal 
transport  barrier  is  in  a minimum  energy  state  under  the  constraint  of  the  conservation  of 
total  angular  momentum. 
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